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I.  INTRODUCTION 


This  is  a  report  of  an  investigation  of  the  feedback  control,  of  a 
variable  air  volume  (VAV)  heating,  ventilating,  and  air  conditioning  (HVAC) 
system  for  a  moderate  size  muLtiroom  building.  The  primary  objective  is  to 
maintain  the  room  temperature  within  a  given  band.  This  report  looks  at  one 
method  of  achieving  this  goal.  This  method  utilizes  proportional  control  of 
the  room  temperatures  and  proportional-integral  control  of  three  other  quanti¬ 
ties.  Each  of  these  controllers  is  isolated  from  the  others  and  attempts  to 
control  its  own  quantity  without  concern  of  the  effect  upon  the  other  quanti¬ 
ties — hence  the  name  decentralized  control.  Of  course,  these  local  loops  are 
expected  to  benefit  each  other  by  maintaining  their  quantity  in  control.  The 
purpose  of  this  investigation  is  to  discern  how  well  these  expectations  are 
met;  that  is,  given  a  good  design  of  the  local  loop  controllers,  what  is  the 
resulting  control  when  these  isolated  subsystems  are  interconnected  into  the 
complete  system? 


II.  THE  SISO  DESIGN  OF  CONTROL 


The  HVAC  system  to  be  studied  is  shown  in  Figure  1.  It  is  caLIed  single 
duct  because  there  is  one  main  duct  providing  cool  supply  air  to  all  the 
zones.  (A  zone  is  a  room  or  a  group  of  rooms  for  which  a  single  temperature 
measurement  is  considered  to  reflect  the  temperature  of  the  entire  space.) 

The  main  duct  contains  a  cooling  coil  which  in  warmer  months  chills  the  supply 
air  to  a  desired  setpoint.  The  zones  in  turn  have  valves  to  decrease  the  vol¬ 
ume  of  cool  supply  air  entering  the  zone.  In  the  winter  months,  the  cooling 
coil  is  shut  off;  the  zone  valves  are  at  a  minimum  flow  setting  and  the  zone 
air  temperature  is  controlled  by  a  heating  coil  which  reheats  the  supply 
air.  Each  zone  has  its  own  reheat  coil  and  air  valve. 

The  VAV  description  means  that  the  air  flow  into  the  zones  can  be  varied 
by  the  air  valves,  or  more  commonly,  Che  VAV  boxes.  In  order  to  avoid 
overpressure  in  the  ducts  and  also  to  allow  energy  savings,  the  speed  of  the 
fan  in  the  supply  duct  can  be  modulated,  resulting  in  variable  air  volume 
flow. 


The  remaining  active  ingredient  of  the  system  is  the  mixing  box  which 
mixes  return  air  from  the  zones  and  outdoor  air.  This  air  is  then  routed  to 
the  face  of  the  cooling  coil  in  the  supply  duct.  The  amounts  of  outdoor  and 
return  air  which  are  mixed  depend  upon  the  desired  mixed  air  temperature. 

All  of  the  controlled  quantities  and  the  means  of  effecting  changes  have 
been  mentioned;  along  with  the  symbols  which  will  be  used,  they  are  now  summa¬ 
rized.  The  mixed  air  temperature,  Tam,  is  effected  by  changing  the  fractions 
of  return  air,  fret>  (at  temperature  Tfet),  and  of  outdoor  air,  l-fret  (flt 
temperature  T0ut),  which  comprise  the  mixed  air.  The  cool  air  temperature  in 
the  supply  duct,  Tcool,  is  affected  by  the  mass  flow  of  water  through  the 
cooling  coil,  m  .  The  gage  pressure  in  the  supply  duct,  Pc,  is  affected  by 
the  fan  motor  speed  or  frequency,  f-n«  The  temperature  of  air  in  a  zone,  Tz£ 
(where  i=l,2,3,4,  i.e.,  there  will  be  f'  r  zones),  is  affected  by:  the  mass 
flow  of  supply  air  into  the  zone,  in  the  summer;  and  the  mass  flow  of 

hot  water  in  the  reheat  coil,  m  . ,  in  the  winter.  Again,  these  quantities 

WZ 1 

are,  first,  the  quantities  which  should  be  controlled  or  regulated  in  a  pre¬ 
scribed  manner  and,  second,  the  quantities  which  can  be  commanded  to  change  in 
an  attempt  to  achieve  control  of  the  first  quantitites;  e.g.,  control  of  Tc00^ 
is  achieved  by  commands  to  m 

wc 

These  input-output  pairs  are  the  basis  for  the  local-loop  subsystems.  Of 
course  they  are  not  the  only  inputs  affecting  the  controlled  quantities;  there 
are  also  the  interconnections.  Although  it  will  be  helpful  when  designing  the 
controllers  to  think  of  the  system  as  local-loop  subsystems,  ultimately  the 
interconnections  must  be  considered.  Figure  2  shows  a  block  diagram  of  the 
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entire  system  decomposed  as  independent  subsystems.  (In  this  figure  the  actu¬ 
ators  have  been  included  in  the  subsystem  blocks,  so  in  place  of  mwc  is  P^, 
the  pneumatic  pressure  to  the  valve  controlling  n>wc  into  the  cooling  coil. 
fret  is  replaced  by  Pad*  the  pneumatic  pressure  to  the  actuator  which  sets 

f_„„.  m  ■  and  m  •  are  replaced,  respectively,  by  P.  -  and  P.  • ,  the 

ret  azi  wzi  ’  dzi  vzi’ 

pneumatic  pressures  to  the  VAV  box  and  water  valve  actuators  which  control 
•  *  . 
m  .  and  m  .  . ) 
a zi  wzi 

It  is  now  possible  to  point  out  the  difference  between  the  single-vari¬ 
able  local-loop  subsystems  and  the  multivariable  system.  All  the  local  loops 
are  single-input  single-output  (SISO)  systems  with  the  inputs  and  outputs 
being  scalars;  e.g.,  for  the  cooling  coil,  the  input  is  P^  and  the  output  is 
Tcool*  further,  the  input-output  relation  can  be  described  by  a  single  trans¬ 
fer  function.  Note  that  there  are  seven  SISO  local  loop  subsystems.  Alterna¬ 
tively,  the  subsystems  can  be  unified  into  a  multi-input  multi-output  (MIMO) 
system  with  the  input  and  output  being  vectors;  the  input,  U,  in  summer  is 

u  -  <p„d  Pvc  f,  pdz,  pdz!  pdzj  pd.»)T  <U-» 

and  in  winter, 

“  *  <P.d  Pvc  PVZ 1  P  VZ  3  Pvz  3  P»„)T  (II-2) 

The  output  is  Y, 


“  *  <p.d  Pvc  £, 


ppp 

dzi  dz2  dz3 

(II-l) 

PPP 

P 

(II-2) 

VZl  VZ  2  vz  3 

vz1*' 

P  T  T  T 

rc  *Z1  AZ2  s3 

T  )T 

Z**' 

(II-3) 

now  described  by  a  7  by  7  matrix. 

whose 

Y  ^Tam  Tcool 


elements  are  transfer  functions.  This  matrix  will  be  calLed  the  system  plant 
transfer  function  matrix.  Note  that  the  seven  on-diagonal  elements  of  this 
matrix  are  identically  the  seven  local  loop  SISO  transfer  functions,  and  the 
42  off-diagonal  elements  describe  the  interconnections  between  the  subsys¬ 
tems.  For  example,  the  second  row  first  column  element  of  the  system  plant 
transfer  function  matrix  describes  the  change  in  Tcoo^  due  to  a  change  in  Pa(j. 

It  is  obvious  that  the  MIMO  approach  is  a  more  complete  description  of 
the  entire  system  than  a  description  by  seven  SISO  subsystems  without  inter¬ 
connections.  Of  course  with  this  increase  in  complexity  of  description  come 
more  questions  when  this  MIMO  approach  is  used  for  design  purposes.  For 
instance,  in  a  typical  SISO  feedback  control  system  there  are  four  elements: 
the  plant,  whose  output  is  the  control  variable;  a  feedback  compensator;  a 
summation  device,  which  compares  the  feedback  signal  to  a  reference  or  set- 
points;  and  a  feedforward  compensator  which  provides  the  input  to  the  plant. 

To  achieve  the  desired  control,  the  design  of  two  elements  is  necessary:  the 
feedback  and  feedforward  compensators.  Alternatively,  for  the  MIMO  descrip¬ 
tion  ,  the  feedback  and  feedforward  compensators  contain  49  elements  each. 

So,  although  the  MIMO  design  may  be  able  to  provide  better  control  of  the 
seven  outputs,  the  difficulty  of  design  is  increased  over  the  SISO  approach  by 
7  to  1. 
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In  order  to  achieve  a  workable  design,  this  investigation  will  use  a  SISO 
approach  to  designing  the  controllers.  Moreover,  the  feedback  compensator 
will  include  no  elements  beyond  what  is  necessary  for  sensing  the  plant  out¬ 
put,  resulting  in  the  need  for  seven  feedforward  compensator  designs. 

The  system  under  consideration  is  a  common  HVAC  system,  and  it  is  desired 
that  the  design  be  implemented  with  commonly  available  components.  For  this 
reason  the  industrial  grade  proportional-integral  (PI)  process  controller  was 
selected  as  the  feedforward  compensators.  The  only  remaining  facet  of  the 
design  is  the  selection  of  the  two  PI  controller  gain  values:  the  propor¬ 
tional  control  gain  Kp  and  the  integral  control  gain  Kj.  These  values  will  be 
obtained  once  the  PI  controller  is  attached  to  the  plant  by  using  the  pro¬ 
cedure  devised  by  J.  G.  Ziegler  and  N.  B.  Nichols1  for  selecting  optimum 
settings  for  automatic  controllers. 

Before  proceeding  to  evaluate  the  controller  design  just  described,  it 
may  be  fruitful  to  describe  what  other  types  of  design  are  possible.  Given 
the  system  plant  transfer  function  matrix  description,  it  is  feasible  to 
provide  precompensation,  which  amounts  to  counteracting  the  interconnection 
effect,  resulting  in  a  new  system  plant  with  all  off-diagonal  elements  being 
zero.  At  that  point  the  simpler  SISO  approach  to  controller  design  would  be 
optimal.  Precompensation  would  entail  selecting  on-  and  off-diagonal  elements 
of  the  feedback  and  feedforward  compensator  matrices.  This  approach,  or  any 
centralized  controller  approach  where  all  seven  plant  outputs  can  be  used  to 
select  appropriate  plant  inputs,  would  increase  in  complexity  at  a  quick  rate; 
however,  given  the  improvement  in  components — e.g.,  electronic-to-pneumatic 
transducers — allowing  more  accurate  modeling  and  the  power  of  Direct  Digital 
Control  for  sensing  and  computation,  the  results  might  be  worth  the  work. 

At  present  the  controller  design  is  SISO  and  our  task  is  to  insure  that 
this  design  provides  good  control.  The  next  step  is  to  develop  our  system 
description  by  modeling  the  plant  and  associated  components. 


III.  MATHEMATICAL  MODELING  OF  THE  SISO  SUBSYSTEMS 


l 
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A  preliminary  description  (Chapter  II)  and  sketch  (Figure  1)  of  the 
system  have  already  been  given.  It  is  now  necessary  to  detail  the  system 
description  to  allow  for  the  use  of  computer  simulation  and  frequency  analysis 
techniques.  The  approach  of  the  next  two  chapters  will  utilize  the  SISO-MIMO 
distinction  made  in  the  previous  chapter;  first,  the  SISO  subsystem  dynamics 
will  be  mathematically  described;  then  in  the  next  chapter  the  MIMO  system 
will  be  described  by  filling  in  the  interconnections  in  the  system  plant 
transfer  function  matrix  and  associated  feedback  and  feedforward  compensator 
matrices . 


There  are  seven  SISO  subsystems:  the  mixing  box,  the  cooling  coil,  the 
fan  and  associated  ducts,  and  the  four  zones.  Because  these  four  types  of 
subsystems  have  many  components  in  common  and  utilize  the  same  feedback 
control  design,  a  standard  form  for  all  the  subsystems  will  be  adopted.  This 
form  is  shown  in  the  block  diagram  of  Figure  3.  The  description  will  now 
proceed  by  selecting  a  block  in  this  diagram,  describing  the  block  for  each 
subsystem,  and  then  selecting  another  block. 


A.  The  Plants 

The  first  plant  to  be  considered  is  the  mixing  box  (see  Figures  A  and 
5).  The  mixing  box  plant  is  the  damper  and  ducts  which  mix  the  return  and 
outdoor  airstreams  to  produce  the  mixed  airstream.  The  mixing  box  will  be 
used  in  one  of  two  ways,  depending  upon  the  return  and  outdoor  air  tempera¬ 
tures.  If  the  outdoor  air  is  hotter  than  the  return  air,  the  damper  is  set 

for  a  maximum  return  air  fraction,  fma*,  and  no  control  of  the  return  air 

ret 

fraction  of  the  mixed  air,  ffet,  is  allowed.  (For  ventilation  purposes  a 
minimum  of  15%  outdoor  air  is  required  in  the  mixed  air,  so,  f™a*  =  0.85.) 
Secondly,  if  the  outdoor  air  is  colder  than  the  return  air,  dynamic  control  of 
the  air  damper  or  fret  is  allowed.  The  relation  between  the  mixed  air  tem¬ 
perature,  Tflm,  and  the  return  air  fraction,  fret»  is  given  by. 


Cm  T 
m  am 


=  -T  ♦  f  T  ♦  ( 1— f  )  T 

am  ret  ret  ret  out 
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where: 


Tout  =  the  outdoor  air  temperature 

Tret  =  the  return  a^r  temperature 

t_  =0.1  sec,  is  the  mixing  box  time  constant 
m 

Note  that  even  during  this  mode,  the  return  air  fraction  is  constrained  to  be 
between  zero  and  85  percent. 


The  cooling  coil  plant  is  considered  to  have  three  inputs:  the  mass  flow 
of  chilled  water  through  the  coil,  m  ;  the  mass  flow  of  air  across  the  coil, 
masup’  an<*  the  temPerature  the  incoming  air,  Tam,  the  mixed  air  temperature 
(see  Figures  6  and  7).  The  output  is  the  leaving  air  temperature  which 
becomes  the  cool  supply  air  temperature,  Tcoop  The  description  of  the 
dynamics  is  taken  in  part  from  a  paper  by  J.  R.  Gartner  and  H.  L.  Harrison2 
which  was  further  developed  by  I.  Gondal3  and  also  in  part  from  experiments 
performed  at  the  U.S.  Army  Construction  Engineering  Research  Laboratory  (USA- 
CERL).  The  Cartner  and  Harrison  derivation  begins  with  an  infinitesimal 
cross-section  of  the  tube  which  is  in  crossflow;  develops  differential  equa¬ 
tions  for  the  primary  fluid  (the  chilled  water),  the  tube,  and  the  secondary 
fluid  (the  cool  air);  and  transforms  these  equations  into  transfer  func¬ 
tions.  For  the  subsystem,  the  important  equation  relates  changes  in  m  to 

•  wc 

changes  in  TcooL,  given  by 


«  T  ,  ♦  «  T  ,  ♦  “  T  . 

2  COOl  L  COOl  0  COOl 


8  m  ♦ 
o  wc 


8  m 
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where  the  <*'s  and  8's  depend  in  a  complicated  manner  upon  the  steady-state 
inlet  water,  inlet  air,  and  outlet  air  temperatures;  the  steady-state  air  and 
water  mass  flows;  and  the  dimensions  and  materials  of  the  coil.  The  effect  of 
changes  in  Tcoo^  due  to  changes  in  Tflm  and  mag  will  be  given  when  the 
interconnections  are  detailed. 


The  second  plant  to  be  considered,  the  fan  subsystem  plant,  includes  the 
fan  motor,  the  fan  cabinet,  and  the  duct  system.  (See  Figures  8  and  9.)  The 
motor  is  a  two  pole-pair,  squirrel-cage  induction  motor,  the  pulley  set  has  a 
one-to-one  ratio,  and  the  fan  cabinet  is  of  the  centrifugal  forward-curved 
blade  type.  The  input  to  the  motor  is  the  frequency  of  the  electric  power 
supply,  and  the  output  is  the  rotational  speed  of  the  fan  shaft,  0^,  in 
RPM's.  As  derived  from  a  model  including  damping  and  the  moment  of  inertia  of 
the  fan  and  from  experimental  observations,  the  input-output  relation  is, 


9  £  ( s  ) 

FTiT  = 
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29.17 
1  +  s  ■  t 
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where  tm  =  0.361  sec.  is  the  fan  motor  time  constant  and  s  is  the  frequency 

domain  variable.  To  relate  the  8-  to  m  and  P  ,  the  fan  curves  must  be 

f  as 

used.  When  operating  as  a  subsystem,  the  fan  is  riding  a  "system"  curve  which 
is  determined  by  the  fraction  of  wide-open-cfm,  fwo,  determined  by  the  ducts 
and  zone  damper  positions.  This  relation  will  be  given  as  an  interconnection, 
and  for  now,  f  will  be  assumed  constant.  By  taking  points  from  the  fan 
curve  and  using  least-squares  curve  fits,  the  following  lines  were  obtained: 


P,  [p  ]  =  A  ♦  B-  Qc  [ RPM ) 
t  a  t 
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where: 


A  =  -1466  ♦  1349  -f 
B  =  2.356-2.169  -f 


where : 


m  [-*-  1  =  A  +  B-Q,[RPM] 
as  sec  fl 


A  =  -  0.1252  ♦  0.661  -f  -  0.6684-f 

wo  wo 

B  =  -  1.417  -10~4  ♦  5.224- 10~3 • f 


( 1 1 1 — 5  ) 


So,  for  changes  in  motor  speed,  60^,  the  resulting  changes  in  mass 

flow,  6m  , 1  and  pressure  across  the  fan,  6Pf,  for  constant  f  are  given  by 

aS  I  WU 

the  constants 


F  =  — -I  ,  =2.356-2.169-f 

l  •  1  6f  =0  wo 

60^  wo 


6m 

— ^L,  “  1.417-  10-3  +  5.224- 10~3  •  f 

.-  6f  =0  wo 

60  j  wo 


( 1 1 1— 6 ) 


(III-7 ) 


Here  it  should  be  noted  that  the  duct  friction  losses  from  the  fan  to  the 
pressure  sensor  are  assumed  to  be  appreciable.  It  was  assumed  that  for  a 
pressure  gain  of  698  pa  across  the  fan,  the  pressure  in  the  duct  was  125  pa  or 


=  0.1786 
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The  last  subsystem  plant  is  that  of  the  zones  (see  Figures  10  and  11). 
This  includes  the  zone  air  and  walls,  and  the  reheat  coil,  which  is  used  only 
in  the  winter.  The  zone  walls  and  contained  air  can  be  described  in  the  heat 
flux  equation. 


C  •  T  .  =  Q  +Q-+Q.  +  L. 
zi  zi  ^w  D  l 


(II 1—9 ) 


lThe  units  of  m  are  kg/sec.  For  an  air  density  of  1.229  kg/m3,  1  kg/sec 
corresponds  to  1724  CFM. 
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where: 


C  •  =  Che  heat  capacity  of  the  air  in  the  zone  (700.  kJ/°C) 

Qw  =  the  heat  flux  through  che  walls  to  the  outdoors 
Qd  =  the  heat  flux  from  the  duct  containing  the  reheat  coil  and 
the  VAV  damper 

Qj  =  the  heat  flux  due  to  infiltration 

=  the  constant  internal  load  in  the  zone  due  to  lights,  equipment, 
and  people. 


The  load  in  zone  i,  ,  is  considered  to  be  a  constant  parameter  but  later  its 
value  will  be  chosen  to  obtain  a  desired  steady-state  zone  temperature.  Qj 
and  Qw  are  described  by 
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where: 

TQut  =  is  the  outside  temperature 

v  =  is  the  conductance  through  the  zone  wails  (0.0644  kw/°C) 

Fj  =  is  the  infiltration  rate  (1/3600  sec-1) 


The  flux  QD  is  given  by 
where: 
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m  .  =  mass  flow  of  air  into  the  room  through  the  duct  (in  the  summer 
flZ  1 

this  is  varied  by  the  VAV  damper,  and  set  to  a  minimum  flow 
position  in  the  winter) 

T^  .  =  temperature  of  the  duct  air  as  it  enters  the  zone 


In  the  summer  when  the  reheat  coil  is  not  operating,  -  ^cool’  *n  w^nter 

the  coil  has  a  hot  water  mass  flow,  m^^,  and  the  duct  temperature  is 
described  by  an  equation  derived  in  the  same  manner  as  the  cooling  coil  (see 
Equation  III-2),  so, 


*  .  T.  ♦  *  .T,  ♦  <*  .T.  =  6  m  .  ♦  0  m 

2zi  dzi  izi  dzi  ozi  dzi  l  wzi  o  wzi 
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where  the  <*'s  and  6's  are  constants  which  depend  upon  the  steady-state 
temperatures  and  mass  flows,  and  upon  the  coil  dimensions  and  materials. 


B.  The  Controllers 


Two  types  of  controllers  are  used  in  this  system.  The  first  is  an  elec¬ 
tronic  PI  controller  which  is  assumed  to  react  so  fast  that  any  dynamics  are 
neglected  and  the  output  is  instantaneous.  As  a  function  of  the  error,  e, 


between  a  temperature  T  and  a  setpoint  temperature  TSP,  i.e.,  e  =  Ts"-T,  the 
controller  output  is  C  , 


(III-14) 


C  ft)  =  K  ■  e(t)  +  Kx  • fe(t)dt  (III-13) 

out  p  I  * 

In  industrial  applications  this  is  sometimes  written 

Cout(t>  =  Kp  [e(t)  *  f  Je<c>dtl  (III-14] 

where: 

Kp  =  the  proportional  gain 
Kf  =  the  integral  gain 

Tj  =  Kp/Kj,  is  the  integral  time  constant 

This  controller  is  used  in  the  mixing  box,  cooling  coil,  and  fan  subsystems. 

The  setpoints  which  are  necessary  to  determine  the  error  term  in  the 

controllers  are  determined  in  various  ways.  The  setpoint  temperature  for  the 

cooling  coil,  T^Pq^,  is  determined  by  the  highest  zone  temperature,  T^g^ 

=  max  [ T  . ] ,  as 
1  zi 


18*3°C  if  Thigh  <24*2°C 
113.4-3.93'T  .  .  if  24.2  <  T.  .  .  <25.6°C 

high  high 

12.8°C  if  T.  .  .  >25.6°C 

high 
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This  variable  setpoint  was  used  to  obtain  the  steady-state  values  but  was  not 

used  in  the  step  tests  of  Chapters  V  and  VII.  The  setpoint  for  the  mixing 

box,  T3p,  is  set  in  coordination  with  TSp  .  ,  as, 
am’  cool 


TSP  =  TSp  -  1 . l°C 
am  cool 


(III-16) 


sp 

Finally,  the  setpoint  for  the  fan  subsystem  controller,  Pc  ,  is  a  constant 
throughout  the  study, 

P3p  =  622  pa  (III— 17 ) 

c  r 

The  zone  subsystem  uses  a  different  controller.  Recall  that  the  zone 
controls  its  temperature  in  two  ways:  1)  by  modulating  the  mass  flow  of  air 
down  the  zone  duct  with  the  VAV  box  and  2)  by  letting  hot  water  flow  in  the 
reheat  coil,  warming  the  supply  air  before  it  enters  the  zone.  This  scheme  of 
control  is  performed  by  the  zero  energy  band  thermostat  (ZET)  (see  Figure 
12).  This  name  is  given  because  there  is  a  band  of  temperatures  (68  to  75SF) 
between  which  it  is  said  that  no  energy  is  consumed.  Actually,  when  in  this 


t ,V. V.VV  *■. 
-  ’ /  V 


band  the  air  flow  is  set  to  a  minimum  which  does  comprise  a  heat  flux,  but  the 
name  is  descriptive  and  common  so  it  will  be  used.  Note  that  because  this 
device  is  pneumatic,  its  dynamics  are  not  negligible.  The  device  has  one 
input,  T2£,  and  two  outputs:  the  pneumatic  pressure  to  the  reheat  coil  water 
valve,  Pvzi>  and  the  pneumatic  pressure  to  the  VAV  box  actuator,  P^^.  The 
ZCT  has  five  modes  of  operation: 


1)  for  Tzi  <  18.3°C  (65°F) 


P.  •  *  Pm^n,  P  • 

dzi  ad  vzi 


(111-18) 


where  pmax  =  89.6  kpa,  Pmjn  =  56.2  kpa 
wv  r  ad 


2)  for  18.3°C  <  T  •  <  20.0 

Z  X 


t  -p  .  =  -p  .  ♦  pmin  ♦  (p1"** 
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(T  20.0) 

Pmin)*  _ 

1 _  *  /  i  a  'i/t  / 


wv  '  (18.3  -  20.0] 


where  tcz  *  1.0  sec,  is  the  time  constant  of  the  ZET  and,  =  62.1  kpa 

3)  for  20.0  <  T  •  <  24.2 

z  1 

p  =  pinin  p  pmin 

dzi  ad  vzi  wv 

4)  for  24.2  <  T  •  <  25.6 

zi 

(T  .  -  24.2) 

t  •  P  .  ■  -P  ♦  Pmxn  ♦  (Pm*x  -  Pm*n)-  - =7-^ 

cz  dzi  dzi  ad  ad  ad  (25.6  -  24.2) 


P  =  Pmin 
vzi  wv 

and  5)  for  T  •  >  25.6 

zi 

_  max  „  _min 

P .  -  =  P  .  ,  P  .  =  P 
dzi  ad  vzi  wv 

This  controller  performs  only  proportional  control.  The  gain  is  set  by  the 
difference  in  the  temperatures  between  which  the  ZET  modulates  the  output 


pressure  and  also  by  difference  of  the  extreme  pressures.11  This  is  different 
than  with  the  PI  controllers  where  the  selection  of  the  gain  settings  is  a 
necessary  and  crucial  procedure,  which  will  be  the  topic  of  the  next  chapter. 


C.  The  Sensors 


In  order  that  the  PI  controller  can  calculate  the  error,  it  must  have  an 
electronic  representation  of  the  control  variable.  This  voltage  is  provided 
by  a  sensing  device  having  a  non-negl igible  dynamic  response.  The  transfer 
function  of  all  the  measurement  operations  has  the  form 


Yu(s)  =  1 

TTsT-  1  +  tjS 
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where  Y°  is  the  Laplace  transform  of  the  measured  variable,  Y  is  the  transform 
of  the  variable  which  is  measured,  and  tj  is  the  sensor  time  constant.  Here  Y 

can  be  Tam*  Tcool’  or  V  Y°  can  be  Toam'  Tocol'  or  Poc»  and  cj  can  be  ^m’ 
tgc,  or  tg£  for  the  mixing  box,  cooling  coil,  and  fan  subsystems,  respec¬ 
tively.  This  defines  the  observed  mixed  air  temperature,  Toam,  where  the 
mixing  box  sensor  time  constant  tsm  *  10.0  secs;  the  observed  cold  deck 
temperature,  Tocoi>  where  the  cold  deck  sensor  time  constant,  t#c  -  10.0  secs; 
and  the  observed  duct  static  pressure  PQC,  where  the  static  pressure  sensor 
time  constant  tg£  =0.1  sec.  Actually,  the  static  pressure  sensor  dynamics 
arise  from  the  sensor  and  a  30.5-m  (100-ft)  tube  necessary  to  reach  the 
duct.  This  length  of  tube  also  introduces  a  transport  lag  or  pure  time  delay 
into  the  fan  subsystem  feedback  due  to  the  finite  speed  of  the  propagation  of 
this  pressure  signal.  This  pure  time  delay  has  a  transfer  function  of 
exp(-s-t^d^)  where  ^^£=0.1  sec  is  the  tube  delay  time. 

The  ZET  is  placed  directly  in  the  zone  so  that  there  is  not  a  separate 
measuring  device  for  the  zone  subsystem. 


D.  The  Actuators 

The  actuators  are  the  final  link  in  the  subsystem,  translating  the  elec¬ 
tronic  signal  of  the  controller  into  an  action  in  the  plant.  For  the  mixing 
box  and  cooling  coil,  this  is  a  two-component  operation:  the  electronic  to 
pneumatic  (E/P)  transducer  produces  a  pressure  signal  from  the  electronic 
controller  output,  and  the  pneumatic  actuator  uses  this  pressure  signal  to 
move  an  air  damper  or  water  valve.  These  two  components  are  described  inde¬ 
pendently.  The  E/P  is  described  by 


riThe  gains  associated  with  the  ZET  were  fixed  throughout  the  course  of  this 
study. 
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cuT  1  ♦  t  -s 
ep 

where: 

P(s)  =  the  transform  of  the  output  pressure:  Pad  for  the  mixing  box; 

P,„  for  the  cooling  coil 

C(s)  =  the  transform  of  the  respective  controller  output 
t._  =  2.4  secs,  the  E/P  time  constant 

ep 

The  E/P  time  constant  was  obtained  from  a  frequency  plot  of  peak-to-peak 
pressure  output  vs.  peak-to-peak  sine  wave  voltage  input.  The  break  frequency 
of  this  corner  plot  was  found  by  extrapolating  a  line  with  -20  dB/decade 
slope,  curve-fit  to  the  highest  frequency  data  points  up  to  the  zero-frequency 
gain  level.  The  inverse  of  the  frequency  at  this  intersection  is  the  time 
constant  of  the  single-pole  transfer  function  model,  Eqn.  III-20. 

Further  confirmation  of  this  model  was  obtained  by  applying  a  step-input 
in  the  voltage  and  observing  the  output  with  a  strip  chart  recorder.  This 
procedure  produced  a  value  of  tg  =  2.0  £  0.5  secs  which  includes  the  value 
obtained  by  the  frequency  plot  method. 

The  pressure  output  of  the  E/P  in  the  mixing  box  subsystem  is  connected 
to  the  actuator  and  damper  combination  in  the  duct  or  mixing  chamber.  The 
relation  between  the  damper  position  and  the  pneumatic  pressure  is  described 
by 


ret(s) 

P  TsT 
ep 


rp1"**  _  pmin; 

v  ad  ad  ’ 


1 

1  t 
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where: 


tam  =  3.2  secs,  the  mixing  box  actuator  time  constant 
Fret(s)  =  the  transform  of  fret,  the  return  air  fraction 
Pep(s)  -  the  transform  of  Pgp 

The  constant  tam  was  obtained  by  the  same  frequency  plot  procedure  as  tgp. 

The  pressure  output  of  the  E/P  in  the  cooling  coil  is  connected  to  the 
chilled  water  valve  actuator.  The  valve  and  actuator  dynamics  are  described 
by 


M  (s)  . 

wc  _  1 _ 

P  (s)  1  ♦  t  s' 

ep  ac 


(pmax  _  pmin) 
rwv  WV 
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where: 

mmax  =  3.15  kg/sec  (50  GPM),“*  the  maximum  mass  flow  of  water 
•  wc 

|  through  the  coil 

M  (s)  =  the  transform  of  m 
wc  wc 

t„,.  =  3.2  secs,  the  water  valve  and  actuator  time  constant, 

ac  ’ 

’  Here  it  has  been  assumed  that  the  dynamics  of  the  valve/actuator  combination 

are  the  same  as  the  damper/actuator  combination,  i.e.,  t  _  =  t  . 

sc  sin 


The  E/P  and  actuator  models  have  been  presented  here  as  separate  and 
isolated  components.  This  is  how  the  frequency  plot  experiments  were 
performed.  However,  when  performing  the  step  input  tests  to  obtain  further 
confirmation,  a  slight  overshoot  of  the  final  damper  position  was  observed 
when  the  E/P  was  given  a  step  input  in  voltage.  Again  the  overshoot  was 
small,  but  perhaps  a  better  model  might  combine  the  E/P  and  actuator  into  a 
single  model.  At  present,  the  transfer  function  between  fret  and  cQut  is 
given  by 
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and  perhaps  a  better  model  would  be 


F  (s) 
ret 
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where  the  polynomial  in  the  denominator  may  have  complex  roots  indicating 
oscillatory  behavior.  For  this  report,  the  first  model  was  used. 


Both  the  E/P  and  the  actuators  suffer  from  the  nonlinear  effect  of 
hysteresis.  This  cannot  be  accounted  for  by  transfer  functions  so  it  cannot 
be  included  in  the  frequency  analysis,  but  it  can  be  included  in  the 
simulation  via  a  subroutine  named  "BACKLASH."  The  hysteresis  problem  was  made 
evident  by  Dolan4  (see  Figures  13  and  14)  for  dampers  and  valves,  and  by  the 
present  author  for  E/P's.  For  the  valve  and  damper  the  hysteresis  is 
characterized  by  having  a  backlash  of  0.5Z  of  full  span;  for  the  E/P's  this 
number  is  1Z.  A  plot  to  demonstrate  this  is  shown  in  Figure  15,  which 
shows  m  with  hysteresis  (MDWC)  vs.  m  without  hysteresis  (PMOWC)  during  a 
simulation  with  standing  oscillations.  Mote  the  difference  between  the  paths 
as  PMDWC  increases  and  as  PMODWC  decreases. 


lllThis  conversion  assumes  a  water  density  of  998.  kg/m3. 


R  FLO\ 


The  zones,  because  they  use  the  pneumatic  ZET,  do  not  require  an  E/P 
transducer.  They  still  require  actuators.  The  transfer  function  between  the 
reheat  coil  water  mass  flow  (m  .)  and  the  pneumatic  pressure  (P,  •)  is  given 

W Z 1  1 

by 
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where: 


■max 


m  .  =  0.126  kg/sec  (corresponding  to  2  GPM),  the  maximum  water 
wzi 

mass  flow  in  zone  reheat  coil 


'az 


=  3.2  secs,  the  zone  actuator  time  constant  with  (taz=tam). 


The  relation  between  the  mass  flow  of  air  into  the  zone  and  the  pneumatic 
pressure  applied  to  the  actuator  can  be  described  in  two  steps:  the  first  is 
dynamic,  the  second  is  numerical.  The  transfer  function  between  the  ZET 
output  damper  pressure,  P<jzi»  and  the  zone  damper  fraction  of  full  open,  f^^, 


is 
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where  f*?in  is  the  minimum  zone  air  damper  fraction  (=0.4).  The  relation 


dz 


between  the  damper  fraction  and  zone  air  mass  flow  is  given  by 


•  -max 
m  =  m 
az  az 
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•max 


where  m‘"“~  is  the  maximum  air  mass  flow  into  the  zone.  This  maximum  is  a 
az 


constant  for  the  subsystem  alone  but  varies  when  the  interconnections  are 
introduced. 


The  actuator  block  of  the  fan  subsystem  is  the  motor  speed  controller. 
This  electronic  unit  converts  a  D.C.  voltage  signal  from  the  PI  controller 
into  a  frequency  signal.  As  the  input  voltage  is  varied  from  0  to  10  volts 
D.C.  the  three-phase  frequency  varies  from  0  to  60  Hz  while  the  output  voltage 
is  varied  to  maintain  constant  volts-per-Hz  ratio.  The  manufacturer  describes 
the  motor  speed  controller  by  stating  that  given  a  step  input  in  voltage,  the 
output  frequency  will  ramp  to  the  new  vaLue.  The  immediate  temptation  is  to 
model  this  as  a  simple  integrator  with  transfer  function,  s  1 .  This  does 
describe  the  output  until  it  reaches  the  new  value  where  the  ramp  is  shut  off, 
indicating  a  transfer  function,  s- 1  (1  -  e-sCi)  where  t;  is  the  time  necessary 
to  reach  the  new  value  given  the  preset  ramp  slope.  For  a  single  known  step 
input,  this  is  useful;  but  for  a  simulation  where  the  inputs  are  not  simple  or 
known  prior  to  modeling,  this  transfer  function  is  not  suitable. 


For  the  simulation,  the  manufacturer's  description  was  used  on  the 
discrete  time  level.  The  simulation  uses  a  Runge-Kutta  integration  technique, 
which,  instead  of  producing  a  continuous  line  vs.  time,  results  in  a  staircase 
output.  The  output  is  constant  over  a  length  of  time  called  the  calculation 
interval,  At,  after  which  the  output  jumps  to  a  new  value.  The  motor  speed 
controller  is  described  by  limiting  the  change  in  output  frequency, 

Af  ^ ,  requested  by  the  input  voltage,  Af°-  If  Af°  is  less  than  m/At  ,  it  is 
allowed,  i.e.. 


min  (Af  ,  m/At)  if  Af  >0 


Af 
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max(Af  ,  -m/At)  if  Af  <0 


Because  this  strategy  is  implemented  on  the  smallest  time-scale  level,  it  can 
accommodate  any  input  voltage  variations  as  well  as  any  other  model.  Also,  it 
does  produce  a  ramp  output  which  stops  discontinuously  at  the  new  value 
selected  by  a  step  input. 


For  the  analysis  techniques,  thi?  model,  (III-28),  is  also  of  no  use  and 
a  less  accurate  model  must  be  used.  For  analysis,  the  relation  between  the 
input  voltage,  CQut,  from  the  controller  and  the  output  frequency,  £  ,  is 
given  by 
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where 


F  ( (s ) 


=  the  transform  of  F, 


CQut ( s  )  =  the  transform  of  CQut 


msc 


=  3.0  secs,  the  motor  speed  controller  time  constant 


The  motor  speed  controller  time  constant  value  was  selected  to  be  appropriate 
to  the  slope  used  which  was  m  =  4  Hz/sec. 


This  concludes  the  SISO  description  of  the  four  subsystem  types.  Next, 
the  MIMO  approach  will  be  used  to  combine  the  subsystems  into  a  whole  system, 
malting  many  subsystem  constant  parameters  into  the  time-varying  inputs  by 
describing  the  interconnections. 


IV.  MATHEMATICAL  MODELING  OF  THE  MIMO  SYSTEM 


Having  described  Che  SISO  local  loops,  ic  is  now  possible  Co  puC  all 
seven  subsysCems  Cogecher  and  form  Che  whole  syscem.  This  will  be  done  in  Cwo 
seeps:  firsc,  Che  incerconneccions  beCween  Che  subsysCems  will  be  presented, 
primarily  from  a  simulaCion  poinc  of  view;  second,  Che  necessary  Cransfer 
funecions  will  be  obCained  from  Che  sysCem  plant  signal  flow  graph  in  order  to 
fill  in  che  elemencs  of  che  plant  transfer  function  matrix. 


A.  Subsystem  Interconnections 

It  is  now  necessary  to  obtain  descriptions  of  how  the  subsystems  affect 
each  other.  Quantities  which  were  assumed  to  be  constant  in  the  isolated 
subsystems  will  not  be  made  variable  as  they  are  affected  by  changes  in  the 
other  subsystems.  The  interconnections  will  be  presented  by  how  they  affect 
each  subsystem. 


The  mixing  box  subsystem  is  affected  by  changes  in  the  zone  air  tempera¬ 
tures  through  the  return  air  temperature.  This  is  shown  by  considering  Tret, 


T 

ret 


! Z  m  . •  T  .  )/(E  m  . ) 

azi  zi  azi 

1=1  1=1 


(IV-1) 


Changes  in  both  T  •  and  m  .  can  affect  T  r .  Also  there  is  a  transport  lag 

time  between  the  time  at  which  Tz-  or  <®az£  changes  and  the  time  at  which  Tret 

at  the  mixing  box  changes.  This  is  due  to  finite  mass  flow  rate  and  duct 

lengths  between  the  zone  and  mixing  box.  This  return  duct  delay  time  is 

tji>  =  61.0/m  ),  where  m  is  the  mass  flow  rate  of  air  in  the  return  duct  in 

dir  as  as 

kg/sec.  This  is  based  on  a  return  air  duct  length  of  27.4  m  (90  ft). 


The  cooling  coil  output,  Tc00p  is  affected  by  two  other  variables,  Tam 

and  m  .  The  derivation  by  Gartner  and  Harrison  of  the  transfer  functions  of 
as 

a  tube  in  crossflow  mentioned  earlier  also  includes  the  relations  between  the 
output  (exiting  secondary  fluid  temperature  Tcoo^),  and  the  inputs  (entering 
secondary  fluid  temperature,  Tam,  and  secondary  fluid  mass  flow,  ®ag)»  The 
first  of  these  transfer  functions  was  found  to  be  unusable  because  it  gave  an 
erroneous  initial  reaction  (although  the  final  values  are  thought  to  be 
correct).  The  second  relation  was  used  and  is  given  by  the  transfer  function 


5  +  5  s  ♦  6  s2 

o  1 _ 2__ 
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where  the  <*'s  are  the  same  as  in  III  —2  and  the  6's  depend  upon  the  same 
quantities  as  the  =  's  and  S's. 

Instead  of  using  the  resuLts  derived  by  Gartner  and  Harrison  for  correla¬ 
ting  changes  in  Tcooi  with  changes  in  Tam,  the  relationship  was  derived 
experimentally.  Using  the  full-scale  HVAC  laboratory  at  USA-CERL  under  the 
direction  of  Richard  Rundus,  an  input  of  T  _  described  by 

T  m(t)  =  t!L  •  t1  "  Exp  C t/a)  ]  (IV-3) 

am  am 

was  imposed  on  the  cooling  coil  operating  at  various  steady-state  flow  condi¬ 
tions.  By  recording  the  TcoQ^  output  data  and  assuming  a  model  of 


a  least  squares  curve-fit  was  made  to  the  predicted  Tc00^  response  of 


(IV-4) 
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This  procedure  resulted  in  the  values  of  G  =  0.4  and  b  =  0.25  sec  1  which 
indicate  a  time  constant  of  the  cooling  coil  for  changes  in  incoming  air  of 
4.0  secs. 


It  should  also  be  mentioned  that  the  temperature  of  the  mixed  air  at  the 

face  of  the  cooling  coil  is  1.1°C  (2°F),  warmer  than  at  the  mixing  box  due  to 

the  heat  generated  by  the  fan  which  is  located  in  the  intervening  duct.  There 

is  also  a  time  delay  associated  with  this  ductwork  which  is  t.,  =  20.0/m 

'  dim  as 

which  is  based  on  a  5.5-m  (18-ft.)  length  of  duct. 


The  next  subsystem  to  be  considered  is  the  fan.  That  there  is  a  connec¬ 
tion  between  the  fan  "system"  curve  and  the  zone  damper  positions  via  the 
fraction  of  wide-open-cfm,  fWQ,  was  mentioned  earlier  (see  Chapter  III — The 
Plants).  This  connection  will  now  be  correlated  with  experimental  data. 
Previous  experiments  on  the  USA-CE&  oVAC  lab  have  produced  an  envelope  on  the 
fan  curve  corresponding  to  various  damper  positions.  For  the  particular  HVAC 
system  under  study — a  single  duct  system,  a  corresponding  f  of  67Z  was 
observed  when  four  dampers  were  wide  open  (fjz£=1.0,  i-1,2,3,4);  when  all  four 
zone  dampers  were  in  the  minimum  position  (f^z£  =  0.4,  i=l,2,3,4),  the  fWQ  was 
seen  to  be  50X.  So,  assuming  that  the  f  moves  linearly  with  fjz£,  the 
relation  becomes 

u 

f  =  0.387  +  0.071  •  1  f  .  ( IV-6 ) 
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The  other  points  on  the  fan  curve  envelope  tend  to  lend  credibility  to  this 
model.  The  delay  between  the  change  in  damper  position  and  change  in  fan 
"system'*  curve  is  assumed  to  be  negligible  (<0.1  sec),  because  the  signal  is 
pressure  wave  propagated  at  the  speed  of  sound. 

Having  given  the  relation  between  f  and  it  is  now  necessary  to 

relate  changes  in  f  to  changes  in  m  and  P  .  This  will  be  done  by  the  same 

wo  as  c 

procedure  which  Led  to  equations  III —6  and  III-7.  From  the  fan  curve  the 
following  two  relations  were  obtained. 


P,[P  ]  =  A  ♦  B  •  f 
f 1  a  wo 


(IV-7) 


where : 


A  =  261.8  -  0.429  •  0f  [ RPM J 


B  =  240.9  -  0.393  •  0f  [RPM] 


mas( kg/sec]  =  A  ♦  B  •  fwQ 
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where: 


A  =  0.1889  -  1.392  10  •  0c[RPM] 


B  =  -0.2747  ♦  5.224  -10  •  0f  [RPM] 


So,  for  a  small  change  in  fwo,  6f  ,  the  changes  in  m  and  Pc  can  be  found  by 
using  the  constants: 


f3  “  i*6.=o  *  24°-9  -  °-J929  • 
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F  =  Tc^1  La  =  "0.2747  ♦  5.224  •  10~  -0f 
u  6f  '60r  _  i 

wo  f=0 


(IV-9) 
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where  0°  is  the  steady-state  operating  point  of  the  motor.  So  using  the 

linearized  fan  model  about  the  steady-state  operating  point,  determined  by 

•o  .  ro 
Qc  and  f  , 
f  wo' 


6P  =  F  (f°  )  •  60  ♦  F  (0p  •  6f 

c  l  wo  f  3  f  wo 


6m  =  F  ( f °  )  •  60r  ♦  F  (0°)  •  6f  „ 
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where  6P  and  6m  are  changes  in  P  and  m  from  their  respective  steady- 
C  d9  ^  s 

state  values,  and  F;  and  F2  are  given  by  III-6  and  1 1 1— 7 . 

The  last  subsystem  is  the  zone.  Here  there  are  two  interconnections 
affecting  the  zone  air  mass  flow  m  ..  and  the  zone  duct  temperature  T2(j£. 
Recall  from  the  SISO  subsystem  zone  actuators  description  (Chapter  III, 
Section  D)  the  equation  III-27, 
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It  is  now  necessary  to  give  meaning  to  the  maximum  zone  air  mass  flow, 

•  max  •  •  •  • 

m  .  This  is  the  maximum  m  .  a  zone  would  be  allowed  to  receive  if  its 
az  azi 

damper  were  moved  to  the  wide-open  position,  f,jz£=1.0.  the  suPP^y  air 

mass  flow,  m  ,  were  for  a  moment  considered  independent  of  damper  position, 

•  max  .  f 

then  m  would  depend  upon  m  and  all  the  zone  damper  positions  as 
a  z  as 
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This  insures  that  when  an  individual  zone  requests  more  air  by  changing  its 
damper  position,  the  change  in  mass  flow  into  that  zone  is  achieved  by  reduc¬ 
ing  the  air  mass  flow  into  the  other  zones  and  not  by  changing  the  total 
supply  air  mass  flow;  note, 
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Of  course, m  does  change  when  the  zone  damper  positions  change,  but  this  is 

due  to  the  fan  'system"  curve  change  and  not  the  manner  in  which  m  is  spLit 
•  •  *  max 

into  the  m  . 's.  Further,  the  m  used  to  calculate  m  is  delayed  from 
.  azi  as  as. 

the  m  in  the  fan  subsystem  by  a  length  of  t,,  =  61.0/m  secs  due  to  the 

as  dls  as 

length  of  supply  duct  between  the  fan  and  the  zones. 


When  the  reheat  coil  in  a  zone  is  operating,  the  change  in  m  .  has 

az  i 

another  effect  which  is  described  in  a  manner  similar  to  the  relation  between 
Tcooi  and  *  •  The  transfer  function  between  T^.  and  m^.is  given  by 
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where  the  * ' s  are  the  same  a9  in  equation  I I I— 12  and  the  6's  depend  upon  the 
same  quantities  as  the  <*'s  and  S's.  Also,  if  the  heating  coil  is  on,  any 
change  in  Tcoo^  will  affect  Tzd£  in  a  manner  similar  to  the  way  Tam  affects 
Tcool'  that  i3»  the  description  is 


T  ..(s)  G 

zdi  =  _ z _ 

Tcool<>>  (1  *  ‘hz** 


(IV-16) 


where  t^z  =  4.0  secs  is  the  same  constant  as  in  IV-4.  The  value  of  Gz,  which 
was  obtained  from  G  in  IV-4  with  the  aid  of  the  Gartner-Harrison  transfer 
functions  for  Tco()^  vs.  Tam  and  Tzd£  vs.  TC0Q^ ,  which  have  reasonable  final 
values,  was  calculated  by  using  the  final-value  theorem  of  Laplace  transforms. 
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This  resulted  in  the  value  Gz  *  0.8.  If  the  heating  coil  is  not  in  use, 

then  any  change  in  Tcoo^  will  produce  an  equal  change  in  Tzd^.  That  is,  Tzdj_ 

=  Tcqo^  when  the  coil's  effect  is  absent.  In  either  case,  coil  or  no  coil, 

the  changes  in  TcoQ^  at  the  cooling  coil  are  delayed  by  time  tdls  = 

61.0/m  secs  before  they  reach  the  zone  duct, 
as 

At  this  point,  all  of  the  necessary  relations  have  been  presented  that 
were  required  to  write  a  simulation  program  of  the  entire  HVAC  system.  This 
was  done  with  the  Advanced  Continuous  Simulation  Language  simulation  package 
which  included  the  pure  time  delays,  the  single  pole  transfer  functions,  the 
quadratic  transfer  functions,  and  the  hysteresis.  This  simulation  will  be 
used  in  the  next  chapters. 

B.  Transfer  Function  Matrix  Description 

In  order  to  use  a  frequency  analysis  technique  to  determine  the  sensitiv¬ 
ity  of  the  seven  subsystems  to  interconnection,  it  i9  necessary  to  derive  a 
matrix  transfer  function  description  of  the  system.  The  MIMO  system  will  now 
be  described  a9  having  vector  inputs  and  outputs  and  by  describing  the 
elements  of  the  transfer  function  matrices.  In  order  that  the  effect  of  the 
interconnections  may  be  studied,  the  MIMO  system  plant  transfer  function 
matrix,  P'  (a  7  by  7  matrix)  will  be  decomposed  into  two  parts  as 
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p*  =  p  *■  AP 


(IV-18) 


where  P  is  Che  nominal  plane  containing  the  disconnected  subsystems  and  AP  is 
the  additive  plant  perturbation  (see  Figure  16).  This  figure  also  shows  the 
controller  transfer  function  matrix,  G,  and  the  feedback  transfer  function 
matrix,  H,  which  also  need  to  be  described. 

The  approach  of  this  section  will  be  to  derive  the  matrix  P'  by  using  a 
signal  flow  graph  of  the  system  plant.  The  matrix  P  will  be  given  without 
need  of  derivation  because  its  diagonal  elements  have  already  been  derived  in 
Chapter  II,  and  its  off-diagonal  elements  are  all  zero  by  the  assumption  that 
the  SISO  subsystems  are  isolated  or  non-interacting.  This  will  then  have 
determined  the  perturbation  matrix  by  AP  =  P"-P.  The  controller  and  feedback 
matrices,  whose  on-diagonal  elements  have  also  been  given  in  Chapter  II  and 
whose  off-diagonal  elements  are  zero  because  the  SISO  subsystems  do  not  have 
interconnections  in  sensing  or  control,  will  then  be  given. 

In  order  that  the  system  plant  matrix,  P',  be  derived,  it  is  necessary  to 
draw  the  signal  flow  graph  of  the  system  plant  (see  Figure  17).  Recall  from 
the  first  chapter  that  the  plant  input  vector  U  in  the  summer  is, 


U  -  <P.d  Pvc  £,  Pd.,  Pd«  pd*=  Pd.,> 


and  in  winter,  it  is 


and  the  output  is  always 
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Because  the  use  of  Pvz^  and  Pdzi  are  made  mutually  exclusive  by  the  ZET,  the 
artifice  of  the  zone  mode  matrix  -  M( i , j ) : i=l ,2 ,3,4; j=l ,2  -  was  introduced 
into  the  system  plant  signal  flow  graph  to  eliminate  the  necessity  of  making 
the  input  vector  have  11  elements.  If  the  zone  i  is  modulating  its  water 
valve,  i.e.,  18.3°C  <  Tzi  <  20.0,  then  M(i,l)  =  1  and  M(i,2)  =  0;  if  the  zone 
i  is  modulating  the  air  damper,  i.e.,  24.2  <  T  •  <  25.6°C,  then  M(i,l)~0, 

M( i ,2 )  =  1;  and  for  all  other  T  j,  M( i , l )=M( i ,2 )=0.  By  including  the  zone 
mode  matrix  in  the  transfer  functions  which  have  Pvz£  or  P^^  as  the  input, 
the  matrix  P*  has  the  dimensions  of  7  by  7. 

The  signal  flow  graph  branches  are  labeled  by  the  transfer  function's  P's 
and  T's.  The  input  and  output  nodes  of  these  branches  are  given  in  Table  1. 


Table  1.  System  PLant  Transfer  Functions 


The  transfer  functions  were  derived  in  the  previous  chapter  and  the  first  part 
of  this  chapter.  They  will  now  be  summarized  here: 
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T  - ( s )  =  •  (T°.  -  T°  )  •  EXP(-s •  t  ,,  ) :  j  =  1,2, 3, 4  (IV-27) 
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Eight  of  the  T'3  depend  upon  the  steady-state  zone  temperature  and  must  be 
specified  in  two  modes: 

a)  if  T° .<  20°C 
zj 
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b)  if  T°.  >  20°C 


m  .  ■  EXP(-s  •  t  ) 

T-  .  (3) - 2^ - - - 4ii - .  j  =  1,2, 3, 4 

11  <C2j'S  *  Vj  * "  *  W 


T  .(•) 

24  +  J 


T°  , .  -  T° 


(C  .  3  ♦  m°  .  +  u  ♦  C  •.  *Ft) 
zj  azj  zj  I 


(IV-33) 


=  1,2, 3, 4  ( IV-34 ) 
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Having  presented  the  detail  of  the  transfer  functions,  it  is  now  possible 
Co  derive,  using  Mason's  gain  formula5,  the  elements  of  the  matrix  P' .  In 
order  to  contain  the  size  of  the  expressions,  a  few  auxiliary  functions  will 
be  defined.  First,  the  loops  on  the  plant  signal  flow  graph  will  be 
designated  by  L's. 
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The  denominator  in  all  of  the  gains  to  be  derived  will  be  D  (sometimes 
labeled  A), 

D(s)  =  1-L1L2LjL4 

The  last  auxiliary  functions  are  the  S'sJ 
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With  the  aid  of  D  and  the  L's  and  S's,  it  is  now  possible  to  present  P',  the 
MIMO  system  plant  transfer  function  matrix  (see  Figure  18). 

The  next  matrix  to  be  presented  is  the  nominal  system  plant  transfer 
function  matrix,  P.  This  matrix,  although  in  MIMO  form,  represents  the  seven 
isolated  SISO  subsystems.  This  matrix  is  shown  in  Figure  19. 
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The  additive  perturbation  matrix  AP  has  now  been  determined.  Using  P' 
and  P, 


AP  =  P'-P 


(IV-42) 


performed  by  simple  matrix  addition. 


The  two  remaining  matrices  are  the  feedforward  and  feedback  transfer 
function  matrices,  G  and  H,  respectively.  These  are  presented  in  Figures  20 
and  21.  The  symbols  in  columns  (or  rows)  four  through  seven  of  H  need  to  be 
defined.  The  AP  are  related  to  the  limits  of  the  pneumatic  pressure  as 
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The  symbol  AT^ 
ZET,  as 


is  associated  with  the  width  of  the  temperature  bands  in  the 


-1.7°C  =  (18.3-20.0)  if  M( j , 1 )=1 
1.4°C  =  (25.6-24.2)  if  M(j,2)=l 


(IV-44) 


This  concludes  the  derivation  of  the  MIMO  description  using  transfer 
function  matrices.  This  description  will  be  used  in  Chapter  VI  to  implement  a 
frequency  domain  analysis  technique  to  determine  the  effect  of  the 
interconnections,  as  a  perturbation,  upon  the  sensitivity  of  the  control  of 
the  system  output. 


C.  Modes  of  Operation!  The  Steady-State  Values 

From  the  presentation  of  the  transfer  function  matrices  it  is  obvious 
that  the  steady-state  values,  e.g.,  T°  _ ,  m°  ,  f°  . ,  are  crucial  to  the 
description  of  the  system.  In  addition,  the  coefficients  of  the  coil  transfer 
functions,  the  «,  0,  and  6's,  depend  upon  the  steady-state  values.  And  when 
the  simulations  are  performed,  Che  steady-state  values  become  so  important 
that  when  a  coil  is  not  in  operation,  it  must  be  deleted  from  the  computer 
input  deck  lest  the  simulation  would  diverge,  producing  erroneous  results. 

In  an  attempt  to  order  the  steady-state  values  in  a  manner  relevant  to 
the  operation  of  the  HVAC  system,  a  mnemonic  was  invented  which  reflects  the 
outdoor  temperature,  Tout,  whether  or  not  the  cooling  coil  is  in  operation, 
and  how  many  of  the  zone  temperatures  fall  in  each  of  the  bands  of  the  ZET. 
This  device  will  be  called  the  "mode  of  operation"  and  will  have  the  form 
mXXLYZ,  where  m  indicates  that  this  is  the  mode  of  operation  mnemonic,  XX  is 
the  outdoor  temperature,  L  is  a  logical  symbol  indicating  whether  the  cooling 
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coil  is  operating  or  not  (it  can  take  values  of  'Y'-Yes,  the  coil  is  operating 

or  'N'-No,  the  cooling  coil  is  not  operating),  and  YZ  which  indicates  how  many 

zones  are  heating-Y  or  cooling-Z  (where  zone  i  heating  means  T°.<20.0  °C  and 

o  z  1 

zone  i  cooling  means  T  ^  >  24.2).  The  sum  of  Y  and  Z  need  not  equal  four;  m 

fact,  the  number  of  zones  in  the  "zero  energy  band"  is  precisely  4-(Y+Z). 

The  selection  of  which  modes  of  operation  were  to  be  studied  was  made  in 
an  attempt  to  represent  those  modes  of  operation  that  an  actual  HVAC  system  of 
this  type  would  require  for  a  climatic  year.  The  modes  of  operation  which 
were  studied  are  presented  in  Table  2. 

There  are  two  subsystems  whose  operational  status  are  not  reflected  in 
the  mode  of  operation  mnemonic:  the  fan  and  the  mixing  box.  The  fan  is 
always  operating,  so  inclusion  of  its  operation  would  be  redundant.  The 
status  of  the  mixing  box  was  left  to  being  more  or  less  a  byproduct  of  the 
mode  of  operation  development  procedure. 

The  modes  of  operation  were  developed  by  selecting  an  outdoor  temperature 

and  the  number  of  zones  which  would  be  heating  or  cooling.  This  determined 

what  the  desired  zone  temperatures  should  be  which,  via  the  cold  deck  reset, 

o 

would  determine  the  steady-state  cold  deck  temperature.  So,  knowing  T  ., 
o  o  ....  . 

Tcqo1»  and  ^out*  and  by  setting  the  derivative  in  III-9  equal  to  zero,  it  was 

possible  to  find  the  internal  loads  of  the  zones,  the  L- 's,  that  were  needed 

to  maintain  T°. .  Knowing  T°.  also  allowed  for  determination  of  T°  .  Because 
zi  ziSD  o  .  •oret 

the  fan  setpoint  was  always  P  -622  pa,  knowing  T  .  determined  m  through 
o  o  ^  o  o  , 

fj  .and  f  .  By  this  procedure  the  resulting  f  ,  .  and  f  were  determined, 
dzi  wo  1  r  dzi.  wo 

Further,  if  the  mixing  box  could  achieve  its  setpoint,  the  cooling  coil,  by 

the  setpoint  selection  procedure,  was  not  required  to  operate.  Because  of 

this,  it  was  not  possible  to  have  the  mixing  box  and  cooling  coil  operating 

simultaneously. 

The  purpose  for  defining  the  mode  of  operation  mnemonic  was  to  improve 
organization.  This  will  be  seen  in  the  following  chapters  where  the  steady- 
state  values  which  were  determined  by  the  whole  system  will  be  used  in  the 
individual  subsystems.  This  insures  that  the  tuning  of  the  subsystems  was  not 
performed  vacuously  but  rather,  with  respect  to  the  values  the  interconnec¬ 
tions  will  have  when  operating  with  the  total  system. 


MODE  tout 


m40y04  40.  25.5  1.0  25.5  1.0  25.5  1.0  25.5  1.0  622  0.662  4.210  13. 


V.  IMPLEMENTING  THE  CONTROL  DESICN:  SELECTION  OP  THE 

CONTROLLER  PARAMETER  VALUES  BY  THE  ZIEGLER-NICHOLS  TECHNIQUE 

Once  the  PI  controllers  have  been  installed  on  the  HVAC  system,  the 
problem  becomes  how  to  determine  the  gain  settings  for  some  "optimum" 
response.  A  practical  method  which  takes  account  of  the  variations  from 
system  to  system  by  performing  a  simple  experiment  is  the  Ziegler-Nichols 
technique. 

J.  Ziegler  and  N.  Nichols1  in  a  1942  paper  entitled  "Optimum  Settings  for 
Automatic  Controllers"  proposed  "Practical  names  and  units  of  measure¬ 
ment... for  the  two  principal  characteristics  affecting  [the]  controllability 
[of  industrial  processes]."  Although  the  authors  speak  of  automatic  reset 
control,  this  is  equivalent  to  the  present  terminology  of  integral  control. 

The  technique  is  performed  by  first  turning  off  the  reset  or  integral 
control  and  then  finding  the  proportional  control  value  entitled  the  "Ultimate 
Sensitivity"  (Su),  above  which  the  output  oscillations  will  continue  to 
increase  in  amplitude  and  below  which  oscillations  will  diminish  to  straight 
line  control.  At  ultimate  sensitivity,  the  amplitude  ratio,  the  relative 
amplitude  of  any  wave  to  that  of  the  wave  which  preceded  it,  is  maintained  at 
unity.  From  the  value  of  Su  and  the  period  of  the  oscillations  (Pu)  at  that 
gain  setting,  the  authors  proposed  that  optimal  settings  for  Kp  and  Kj  can  be 
derived.  The  optimality  of  the  settings  is  a  qualitative  tradeoff  resulting 
in  "reasonably  rapid  [response]  without  excessive  loss  of  stability  or 
excessive  increase  in  period." 

The  steps  to  implement  this  technique  are  clearly  presented  in  the 
literature  which  is  provided  by  the  controller  manufacturers. 

1.  Lock  out  integral  and  derivative  control. 

2.  Set  proportional  gain  Kp  at  sufficiently  small  value  to  rule  out  any 
oscillation. 

3.  Increase  Kp  until  the  control  variable  starts  to  oscillate,  i.e., 
oscillation  once  started  does  not  die  down  but  is  sustained  (this  value  of  Kp 
is  Su). 


4.  Measure  the  period  of  oscillation  at  this  setting  (this  is  Pu)  (see 
Figure  22). 
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5.  The  optimum  values 
the  following  formulas: 


for  K„  and  Kr 
P  1 


can  then  be  calculated  according  to 


K  =  0.45  •  Su  (V-l) 

P 

K  =  0.54  •  Su/Pu 


or  if  the  integral  time  constant  is  used  it  is 

Tj  =  0.83  Pu  (V-2) 

This  technique  for  control  parameter  selection  will  be  used  on  the  mixing 
box,  cooling  coil,  and  fan  subsystems.  The  steady-state  data  will  be  used  for 
the  constant  input  values  which  would  be  variable  if  the  subsystems  were 
interconnected.  The  effects  of  seasonal  changes  of  the  subsystem  components 
on  the  optimum  parameter  values  can  then  be  seen.  Also,  the  effect  of 
hysteresis  will  be  shown  by  finding  optimum  parameter  values  for  both  the 
linear  and  non-linear  subsystems.  Finally,  it  will  be  found  that  Equation 
(V-l)  results  in  a  response  which  is  highly  oscillatory  or  underdamped, 
whereas  a  critically  damped  or  1/4-wave  response  is  more  desirable.  By  chang¬ 
ing  the  fractions  in  (V-l),  parameter  values  resulting  in  this  l/4-wave 
response  can  be  found. 

The  first  subsystem  to  be  presented  is  the  mixing  box.  The  control 
action  is  operative  during  three  of  the  modes  of  operation:  mlOnOO,  ml0n21, 
ml5n04.  The  results  of  performing  the  Z-N  test  on  these  three  sets  of  steady- 
state  values  are  shown  in  Table  3.  From  these  results  it  can  be  seen  that  the 
introduction  of  an  E/P  and  a  damper  with  hysteresis  reduces  the  ultimate 
sensitivity  by  26Z  and  slightly  more  than  doubles  the  ultimate  period,  Pu. 

The  change  in  mode  of  operation  changes  Su  at  most  by  15%  and  the  effect  on  Pu 
is  less  than  6%.  More  dramatic  is  the  need  to  reduce  the  gains  to  achieve  a 
1/4-wave  response.  In  the  linear  case,  the  reduction  is  80%  and  in  the 
nonlinear,  60%.  Plots  of  the  response  of  TQam  to  a  1°C  increase  in  the 
setpoint  for  the  linear  and  nonlinear  systems  with  both  underdamped  and  criti¬ 
cally  damped  control  parameter  settings  for  the  mlOnOO  mode  of  operation  are 
shown  in  Figures  23,  24,  25,  and  26.  From  the  results  of  these  tests  a 
revised  formula  for  obtaining  the  settings  necessary  for  a  1/4-wave  response 
from  the  linear  system  is 


K  =  0.09-S  (V-3) 

P  u 

Kt  =  0.11-S  /P 
I  u  u 
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HYSTERESIS 


TIME  (sec) 


Nonlinear  mixing  box  secpoint  step  change 
-  critically  damped  response. 


For  the  nonLinear  system,  a  quarter-wave  response  could  be  obtained  by  using 


K  =  0.18-S  (V-4) 

P  u 

Kt  =  0.22-S  /P 

I  u  u 

The  next  subsystem  is  the  cooling  coil.  The  coil  i9  in  operation  during 
four  modes  of  operation:  m20y00,  m20y03,  m30y04,  m40y04.  The  results  of 
performing  the  Z-N  test  on  these  sets  of  steady-state  values  are  shown  in 
Table  4.  From  these  results  it  can  be  seen  that  the  introduction  of 
hysteresis  in  the  E/P  and  the  water  valve  reduces  Su  by  10,  16,  17,  and  18% 
and  increases  by  59,  41,  40,  and  27%  for  the  modes  m20y00,  m20y03,  m30y04, 
and  m40y04,  respectively.  The  changes  in  Su  and  Pu  for  different  modes  of 
operation  are  quite  large;  the  maximum  change  for  Su  is  more  than  200%  and  for 
Pu  more  than  20%.  Again,  there  was  a  need  to  drastically  reduce  the  control 
parameter  settings  from  the  value  suggested  by  Z-N  to  achieve  a  1/4-wave 
response.  A  60%  reduction  was  necessary  for  the  linear  system.  A  50% 
decrease  was  needed  for  the  system  with  hysteresis.  This  suggests  a  revised 
formula  set  to  obtain  a  critically  damped  response  to  a  step  change  in  the 
setpoint.  For  the  linear  system,  they  are: 

K  =  0.18-S  (V-5) 

P  “ 

Kt  =  0.22-S  /P 

I  U  U  ' 

and  for  the  non) inear  system, 

K  =  0.23-S  (V-6) 

P  u 

Kt  =  0.27-S  / P 

I  u  u 

Plots  of  cooling  coil  output,  TQC0^,  for  a  step  change  of  1°C  in  the  setpoint 
when  the  system  is  operating  in  the  mode  m20y03  for  the  Linear  and  nonlinear 
components  and  underdamped  and  critically  damped  settings  are  shown  in  Figures 
27,  28,  29,  and  30. 

The  next  subsystem  is  the  fan.  The  fan  is  in  operation  during  all  the 
modes,  although  four  modes  have  the  same  steady-state  conditions  for  the 
fan.  The  results  of  performing  the  Z-N  test  on  the  different  sets  of  steady- 

state  values  are  shown  in  Table  5.  For  this  subsystem,  the  maximum  change, 

due  to  different  modes,  is  14%  for  S  whereas  Py  is  the  same  for  all  the 
modes.  The  need  to  reduce  the  gains  to  achieve  a  1/4-wave  response  was  again 
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Figure  27, 
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Linear  cooLing  coil  setpoint  step  change 
-  underdamped  response. 
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Linear  cooling  coil  setpoint  step  change 
-  critically  damped  response. 


present.  Here,  the  reduction  was  80Z.  So,  the  revised  formula  for  achieving 
a  1/4-wave  response  in  the  isolated  fan  subsystem  is 

K  =  0.09-S  (V-7) 

P  u 

Kt  =  0.11-S  / P 
I  u  u 

Plots  of  the  response  of  PQC  and  mflg  to  a  step  change  in  setpoint  for  the  mode 
of  operation  m30y04  are  shown  in  Figures  31  and  32  for  the  underdamped  case 
and  in  Figures  33  and  34  for  the  critically  damped  case. 

The  final  subsystem  is  the  zone.  This  subsystem  is  controlled  by  the 
zero  energy  band  thermostat  (ZET)  which  is  a  proportional  controller.  The 
gains  are  set  by  the  widths  of  the  temperature  bands  in  which  the  valve  and 
damper  are  modulated.  The  temperatures  at  the  extremes  of  these  bands  are  set 
at  the  factory.  Further,  because  these  controllers  are  relatively  cheap  and 
serve  an  energy  conservation  purpose  beyond  that  of  dynamic  control,  they  were 
the  only  controllers  for  the  zones  considered  in  this  study. 
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Figure  33.  Fan  critically  damped  response  to  a  setpoint  step  change 
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VI.  THE  EFFECTS  OF  INTERCONNECTION:  SINGULAR  VALUES 

BOUNDS  FOR  COMPARISON  SENSITIVITY 

This  chapter  will  present  a  technique  for  determining  whether  a  condition 
on  the  sensitivity  of  the  SISO  design  is  met  when  the  interconnections  are 
considered  to  be  present.  The  condition  on  the  sensitivity  is  called 
comparison  sensitivity.  When  the  comparison  sensitivity  condition  for  the 
system  is  met,  the  sensitivity  to  modeling  errors  of  the  closed  loop  system  is 
less  than  for  an  open-loop  system  which,  without  the  modeling  errors,  has  the 
same  overall  transfer  function  as  the  closed-loop  system.  So,  if  the 
comparison  sensitivity  is  not  met,  there  exists  an  open-loop  system  with  the 
desired  (by  design)  overall  transfer  function  which  is  less  sensitive  to 
modeling  errors  than  the  designed  closed-loop  transfer  function  of  the  total 
system.  Obviously,  comparison  sensitivity  is  a  desired  goal  of  a  closed  loop 
system.  In  fact,  Bode,  a  pioneer  in  the  design  of  feedback,  systems, 
considered  compensation  for  unmodeled  dynamics  to  be  the  primary  motivation 
for  closed-loop  control.6  For  a  more  technical  derivation  of  the  comparison 
sensitivity,  see  Cruz,  Freudenberg,  and  Looze.7  Leave  it  to  be  said  that  for 
comparison  sensitivity  to  be  met  is  to  assure  that  the  inequality 

f  ”  e^(t)e  (t)dt  <  J  ”  e  (t)e  (t)dt  (VI-1) 

J  O  C  C  J  O  0  o 

holds,  where  ec(t)  is  the  error  vector  of  the  difference  between  the  modeled 
plant  and  the  true  plant  for  the  closed-loop  system  and  eQ(t)  is  the  error 
vector  for  the  open-loop  system. 

The  errors  in  modeling  a  plant  are  commonly  described  in  two  ways: 
multiplicative  or  additive  perturbations  to  the  nominal  plant  model.  For  the 
present  system,  the  additive  perturbation  description  will  be  used  (see  Figure 
16).  It  is  now  necessary  to  have  a  test  for  comparison  sensitivity  when  the 
plant  perturbations  are  additive.  This  is  given  by  Cruz,  et  al.9,  using 
singular  values.1  The  test  is: 

Theorem:  Consider  a  nominal  plant  P,  the  additive  plant  perturbation  AP,  the 
feedback  compensator  H  and  the  feedforward  compensator  G.  Assume  that 
(I+PGH)-1  and  AP  are  asymptotically  stable.  If  AP  satisfies 


1The  notation  for  a  singular  value  of  the  matrix  A  is  given  by  cm  [A] 
is  the  positive  square  roots  of  the  eigenvalues  of  AHA  (where  AH  ; “ 
the  conjugate  transpose  of  A). 


is 


and 


min  ck  ( I  +  P(  jui)G(  jw)H(  jw)  ] 

>  max  o.  [  AP(  jto)  ] -max  o.  [G(  ju>)H(  ju>)  ]  +  1  (VI-2) 

for  all  and  the  additiveLy  perturbed  system  is  asymptotically  stable 

for  u)>  a)^ ,  then  the  additively  perturbed  system  is  asymptotically  stable  and 
satisfies  the  sensitivity  condition  (VI-1)  for  all  reference  inputs  r(t)  such 
that 

R(ju)  =  0  for  uj  >  a)  (VI-3) 

c 

where  R(s)  is  the  transform  of  r(t). 

An  interpretation  of  VI-2  is  that  the  smallest  singular  value  of  I  +  PGH 
-  the  distance  in  n-space  to  the  nearest  singular  matrix,  where  n  is  the 
largest  dimension  of  the  matrices  -  must  be  larger  than  the  product  of  the 
largest  singular  values  of  AP  and  GH  -  the  distance  to  the  furthest  singular 
matrices  -  plus  an  amount  necessary  to  insure  the  sensitivity  condition. 

(Note  that  without  the  addition  of  unity,  VI-2  would  express  the  condition 
necessary  to  insure  asymptotic  stability  of  the  perturbed  system.)  It  is 
important  to  note  that  the  test  is  sufficient  but  not  necessary.  If  the 
inequality  in  VI-2  is  satisfied,  the  comparison  sensitivity  is  guaranteed;  the 
converse,  however,  does  not  hold.  If  the  test  fails  (the  inequality  in  IV-2 
is  not  met  for  some  uj  <  ui^),  the  comparison  sensitivity  may,  in  fact,  hold. 

Now  that  the  test  has  been  presented,  it  will  be  used  on  the  system  under 
study.  Recall  that  the  matrices,  P',P,  AP,  G,  and  H  have  already  been  derived 
(see  Chapter  IV).  Using  the  steady-state  data  of  the  various  modes  of  opera¬ 
tion  and  the  controller  parameter  values  of  Chapter  V,  the  quantities  on  the 
right-  and  left-hand  sides  of  VI-2  must  be  computed  and  compared. 

For  some  modes  of  operation,  some  of  the  columns  (or  rows)  of  G  will 
contain  all  zeros.  This  happens  because  one  of  the  control  variables  is 
uncontrolled  for  Che  specific  mode  of  operation,  e.g.,  when  Tz^  falls  in  the 
dead  band,  the  cooling  coil  is  not  on,  or  TQut  >  T  which  sets  fret  = 

0.85.  Because  this  results  in  both  sides  of  VI-2  equaling  unity,  the  test  is 
of  no  value.  In  order  to  use  the  test,  only  the  variables  which,  for  the 
present  mode,  are  being  controlled  should  be  considered  as  elements  of  the 
output  vector  Y.  This  is  done  by  deleting  rows  and  columns  from  P,  P' ,  AP,  G, 
and  H.  Note  that  because  all  of  the  forward  path  gains  were  considered  in 
deriving  P',  deletion  of  the  row  and  column  associated  with  the  mixing  box 
when  it  is  not  controlling  will  not  result  in  deletion  of  the  effect  of 
changes  in  TreC  on  T  p  Likewise,  this  is  true  of  deleting  any  row  and 
column  pair  when  it  is  not  controlling. 
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So,  for  Che  mlOnOO  mode,  where  the  cooling  coil  and  the  zones  are  not 
controlling,  the  matrices  have  a  dimension  of  2x2  corresponding  to  the  fan  and 
mixing  box  working.  The  results  of  the  comparison  sensitivity  test  are  shown 
in  Figure  35.  For  ml5n04  the  matrices  are  6x6  because  the  cooling  coil  is 
off.  Results  of  the  test  for  this  mode  are  shown  in  Figure  36.  For  m20y00 
the  matrices  are  2x2  and  the  test  results  are  shown  in  Figure  37.  The  M40y04 
mode  has  6x6  matrices  and  the  test  results  are  shown  in  Figure  38. 

By  inspecting  the  Figures  35,  36,  37,  and  38  and  referring  to  Equation 
VI-2  it  can  be  seen  that  the  test  failed  on  every  try.  In  fact,  the  system  in 
every  mode  of  operation  failed  to  satisfy  this  test  for  comparison  sensitiv¬ 
ity.  Of  course,  this  does  not  mean  that  the  system  necessarily  fails  to  meet 
the  comparison  sensitivity  condition  (Equation  VI-1);  this  test  simply  fails 
to  prove  it. 

These  results  should  not  be  construed  to  imply  that  the  test  is  not 
useful.  In  fact,  an  option  is  to  now  consider  this  test  as  a  design  criter¬ 
ion,  i.e.,  devise  new  G  and  H's,  possibly  with  non-zero  off-diagonal  elements, 
which  along  with  P  and  AP,  satisfy  the  test  given  in  Equation  VI-1.  The 
procedure  would  then  insure  that  the  comparison  sensitivity  condition  was  met. 

One  method  of  proceeding  to  find  new  G  and  H's  would  be  to  keep  the  SISO 
design  of  the  controllers  in  G  but  also  include  precompensation  in  G  and  H 
which  would  anticipate  the  perturbative  effects  of  AP,  the  interconnections. 
Although  this  may  '•>  analytically  possible  and,  with  the  advent  of  Direct 
Digital  Control  in  HVAC  systems,  be  technically  feasible,  it  may  lead  to  a 
highly  conservative  design.  This  has  been  a  criticism  of  singular  value 
techniques  in  the  past.9 

Before  we  embark  upon  this  formidable  task,  perhaps  the  present  SISO 
design  can,  by  another  method,  be  shown  to  result  in  good  control. 


v  v  v  v  v  ' 


%’VVVV'  •/ 


ran 


SINGULAR  VALUE  PLOT 


Comparison  sensitivity  test  -  mlOnOO 
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sensitivity  test  -  m40y04 


VII.  THE  EFFECTS  OF  INTERCONNECTION:  COMPUTER  SIMULATION 


This  chapter,  like  the  previous,  will  attempt  to  discern  the  effects  of 
the  presence  of  the  interconnections  on  the  SISO  designed  control  of  the 
subsystem  controL  variables.  Having  seen  the  comparison  sensitivity  technique 
fail  to  establish  knowledge  of  these  effects  for  our  present  system,  it  is 
necessary  to  choose  a  different  approach.  This  approach  will  utilize  the 
computer  simulation  developed  in  Chapters  III  and  IV  along  with  steady-state 
modes  of  operation  and  the  controller  parameter  values  found  in  Chapter  V.  As 
in  Chapter  V,  step  changes  will  be  made  in  the  setpoints  of  the  mixing  box, 
cooling  coil,  and  fan  subsystems;  additionally,  load  changes  will  be  made  to 
occur  in  the  zones.  By  recording  the  response  of  the  control  variable  in  the 
primary  subsystem1  and  the  variations  of  the  secondary  subsystems11  control 
variables,  due  to  the  interconnections  to  the  primary  subsystem,  it  will  be 
possibLe  to  see  the  effects  of  interconnection,  as  calculated  by  the 
simulation.  The  effect  of  the  interconnection  will  be  seen  in  two  ways. 

First,  by  comparing  the  response  of  the  primary  subsystem  to  the  isolated 
subsystem  response,  it  will  be  possible  to  observe  the  degradation  of  the 
response  from  the  1/4-wave  response.  Second,  any  variation  in  the  secondary 
subsystem  control  variables  will  be  different  from  the  SISO  design  where,  by 
assuming  the  interconnection  variables  to  be  constant,  it  was  assumed  that  the 
secondary  subsystem  control  variables  maintained  their  constant  steady-state 
values.  Again,  the  SISO  design  of  the  controllers  will  be  tested  for  MIMO 
quality  by  observing:  1)  the  degradation  of  the  primary  subsystem  response  to 
a  step  change  and  2)  the  variations  of  the  secondary  subsystems  control 
variables . 

The  chapter  will  proceed  by  considering  step  changes  in  the  setpoints  of 
the  mixing  box,  the  cooling  coil,  and  the  fan.  Then  a  step  change  will  be 
made  to  occur  in  one  of  the  zones.111 

A.  Mixing  Box  Setpoint  Step  Change 

The  mixing  box  underwent  a  setpoint  step  change  during  two  modes  of 
operation:  mlOnOO  and  ml5n04. 

A  summary  of  the  result  of  these  simulations  is  shown  in  Table  6.  For 

s  p 

ml5n04,  the  step  change  to  the  setpoint  was  ATfl^  =  0.5;  there  was  no  degrada¬ 
tion  from  design  of  the  primary  subsystem  response:  the  zone  temperatures 


1Primary  subsystem:  The  subsystem  which  incurs  the  step  change  in  set 
, , point  or  load. 

.  “Secondary  subsystem:  A  non-primary  subsystem. 

111The  reset  of  the  cooling  coil  setpoint  by  the  highest  zone  temperature  was 
not  considered  in  this  chapter. 
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Mixing  Box  Setpoint 


(all  responding  identically)  changed  by  0.32,  the  pressure  in  the  duct  varied 
<0.12lv,  and  the  supply  air  mass  flow  increased  by  42.  Overall,  for  this 
mode,  the  variations  were  very  small  and  there  was  no  degradation.  (It  should 
be  noted  that  mflg  is  not  a  control  variable  but  is  presented  because  it  is  the 
link  between  the  fan  and  the  other  subsystems.) 

For  the  mlOnOO  mode,  the  only  secondary  subsystems  affected  were  the 
zones.  This  was  because  the  zones  were  in  the  zero-energy  band  and  were  not 
modulating  their  dampers.  The  change  in  T  .  was  -22,  but  from  energy 
considerations,  this  is  considered  acceptable.  There  was  a  slight  degradation 
of  the  mixing  box  step  response.  To  see  this,  observe  Figure  39,  for  the 
isolated  subsystem  response;  Figure  40,  for  the  primary  subsystem  response; 
and  Figure  41,  for  the  secondary  subsystem  response. 

This  concludes  the  mixing  box  setpoint  step  change  simulations.  In 
summary,  a  change  in  the  mixing  box  subsystem  seems  to  have  either  minimal  or 
acceptable  effects  on  the  secondary  subsystems. 

B.  Cooling  Coil  Setpoint  Step  Change 

The  same  procedure  just  described  for  the  mixing  box  subsystem  was 

performed  on  the  cooling  coil.  The  results  of  the  simulations  are  summarized 

in  Table  7.  For  all  the  simulations,  the  step  change  had  a  magnitude  of 

unity.  Most  of  the  primary  responses  suffered  degradation,  due  to  Am  *  0 

because  the  zone  dampers  moved.  The  variance  of  Pc  was  less  than  a  single 

Pascal,  the  maximum  Am  was  -0.200  for  m40y04  with  ATsj*  .  =  -1,  and  the 

as  sec  _  7  3p  cool 

maximum  At  ^  was  +0.4°C  change  for  m40y04  with  ATcj^  =  -1. 

An  example  of  the  effect  of  the  interconnections  between  the  cooling 
coil,  the  zone,  the  fan,  and  back  to  the  cooling  coil  by  virtue  of  TCOQ^  being 
affected  by  a  Am  *0  was  observed  in  the  m20y03  mode  simulation.  A  plot  of 

AS 

the  isolated  subsystem  response  is  shown  in  Figure  42.  The  primary  subsystem 
response  is  shown  in  Figure  43.  The  secondary  subsystem  responses  are  shown 
in  Figures  44,  45,  and  46.  To  obtain  a  measure  of  the  degradation,  the  settl¬ 
ing  time,  tget  (the  time  required  for  the  control  variable  to  stay  within  10 
percent  of  the  final  value),  for  the  two  responses  was  estimated.  For  the 
isolated  subsystem  tget*70  sec3  and  *or  the  primary  subsystem,  secs. 

Note,  however,  that  they  both  eventually  assume  the  setpoint.  Also,  the 
variations,  although  noticeable,  do  not  appear  to  be  significant. 


lvAs  a  point  of  clari f icat ion,  the  S's  in  the  associated  tables  are  the 
maximum  changes  of  the  value  divided  by  its  steady-state  value,  whereas 

the  A's,  e.g.  A  ^  .  are  the  differences  between  the  final  value  and  the 

.  #  .  mas 

initial  value  of  the  variable. 
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Overall,  m20y03,  just  observed,  was  one  of  the  worst  cases,  using 
degradation  and  variation  as  measures,  and  this  mode  was  seen  not  to  suffer 
intolerably  from  the  disruption  of  a  cooling  coil  setpoint  step  change.  It 
seems  fair  to  conclude  that  the  system  is  fairly  insensitive  to  this  change. 

C.  Fan  Setpoint  Step  Change 

The  procedure  used  in  the  previous  two  sections  will  be  used  again  here 
to  examine  the  results  of  a  step  change  of  the  setpoint  for  the  fan 
subsystem.  This  test  was  performed  using  five  of  the  steady-state  modes  of 
operation.  The  results  of  these  tests  are  summarized  in  Table  8.  All  of  the 

S  D 

step  changes  were  the  same,  AP  K  =  50  pa.  The  maximum  of  the  recorded  mass 

c 

flow  changes,  22Z,  occurred  when  the  dampers  in  the  zones  were  not  operating 
(mlOnOO).  The  maximum  change  in  Tz£,  6 Z,  also  occurred  when  fwo  did  not 
change  (m00n40).  The  mixing  box  was  not  operating  for  three  modes,  suffered 
no  change  for  another,  and  varied  6Z  for  mlOnOO,  but  returned  to  the  setpoint 
with  a  settling  time  (10Z  of  maximum  deviation)  of  300  secs.  The  cooling  coil 
was  operating  during  two  of  the  modes  tested,  and,  although  the  cold  deck 
temperature  did  return  to  the  setpoint,  the  variations  were  rather  large,  40Z 
and  43Z.  This  effect  was  due  to  the  change  in  air  mass  flow  across  the  coil. 

As  an  example  of  this  setpoint  change,  the  simulated  outputs  for  m40y04 
are  presented  in  Figure  47  for  P  and  Figure  48  for  m  in  the  isolated 
subsystem.  Figures  49  and  50  for  P  and  m  in  the  primary  subsystem,  and 

v  AS 

Figures  51  and  52  for  the  secondary  subsystems.  The  Pc  responses  were 

indistinguishable;  this  was  the  case  for  every  mode.  The  m  's  have  the  same 

A  9 

initial  response,  but  the  primary  subsystem  does  not  stay  at  the  steady-state 

value  of  isolated  response;  it  again  begins  to  decrease  due  to  the  closing  of 

the  zone  dampers.  The  cooling  coil  is  again  affected  by  Am  but  returns  to 

as 

the  setpoint.  The  zone  temperature  is  lowered  by  the  increase  in  Am  and  the 

A  8 

decrease  in  Tcool. 

In  this  subsystem  setpoint  step  change,  there  is  evidence  of  variations 
in  the  secondary  control  variables.  But,  because  the  control  variables  with 
integral  control  return  to  their  setpoint  and  because  those  with  proportional 
have  only  small  changes,  the  variations  may  be  considered  acceptable.  And 
again,  there  was  no  degradation  of  the  primary  control  variable  response. 

D.  Zone  Load  Step  Change 

To  test  the  effect  of  interconnection  by  changes  in  the  zone  subsystem, 
the  zone  load  L-  was  given  a  step  change  and  the  system  control  variables  were 
observed.  This  test  was  performed  using  seven  of  the  modes  of  operation.  The 
results  of  these  tests  are  summarized  in  Table  9.  The  maximum  change  in  Tz^ 
for  the  smallest  load  change,  ALOAD,  was  9Z  for  mlOnOO.  The  secondary  zone 


Table  8.  Fan  Setpoint  Step  Change  Response 
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Figure  48.  Fan  (m40y04)  setpoint  step  change 
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subsystem  temperature  changes  were  barely  perceptible,  the  maximum  being  1Z 
for  m20y00. 

For  two  of  the  modes,  the  mixing  box  was  operating.  For  mlOnOO  the  Tam 

deviated  by  <  0.9Z  and  returned  to  the  setpoint.  For  ml5n04  the  deviation  was 

so  small  that  the  value  of  the  integral  in  the  controller  did  not  become  large 

enough  to  overcome  the  hysteresis  within  the  duration  of  the  simulation  (1000 

secs);  hence  the  deviation  6  equals  the  changes  AT  in  percentage. 

z 

The  cooling  coil  was  operating  during  four  modes.  For  m20y03,  with  a 

load  change  in  zone  1,  there  was  no  change.  For  m20y03,  with  a  load  change  in 

zone  2,  and  for  m40y04,  the  cold  deck  temperature  deviated  from  the  setpoint 

by  a  small  amount  due  mainly  to  Am  ,  and  then  returned  to  the  setpoint.  For 

^  3 

m20y00,  the  steady-state  chilled  water  mass  flow  is  low,  resulting  in  large 
gains  in  TC0Q^  vs.  Amag  (as  simulated  by  the  transfer  functions).  This  change 
is  so  large  that  Tc00^  could  not  be  forced  to  return  to  the  setpoint  by 
opening  the  water  valve  to  full  flow. 

The  fan  was  operating  in  all  of  the  modes  tested.  For  three  modes  the 
zone  damper  positions  did  not  change  so  the  fan  speed  did  not  change.  For  the 
modes  where  the  fan  was  active,  the  pressure  in  the  duct  deviated  from  the  set 
point  by  <  0.1Z.  The  air  mass  flow  had  large  changes  of  up  to  8.8Z,  but  mag  is 
not  a  control  variable  so  this  change  is  important  only  in  that  it  affects  the 
other  subsystems. 

Two  examples  of  this  test  will  be  presented.  The  first  is  the  mlOnOO 
mode  with  a  -3.0  kW  step  change.  The  plots  are  shown  in  Figure  53  for  the 
isolated  zone,  Figure  54  for  the  primary  zone,  and  Figures  55  and  56  for  the 
secondary  zones.  Comparing  Figures  53  and  54  reveals  that  the  intercon¬ 
nections  had  little  effect  on  the  zone  response.  The  secondary  zone  is  seen 
to  decrease  slightly  due  to  a  decrease  in  Tcoo^  due  to  the  decrease  in  Tz£. 

The  mixing  box  response  is  again  small,  but  the  effects  of  hysteresis  are 
clearly  evident. 

The  second  mode  to  be  presented  is  m40y04  with  a  +3.0  kW  step  change  in 
the  load  of  zone  1.  The  plots  are  shown  in  Figure  57  for  the  isolated  zone. 
Figure  58  for  the  primary  zone  and  Figures  59  and  60  for  the  secondary  subsys¬ 
tems.  Again  the  difference  in  response  for  the  primary  subsystem  is  indistin¬ 
guishable  from  that  for  the  isolated  subsystem.  The  deviations  in  TcOQ^  are 
clearly  evident  and  so  too  is  the  return  to  the  setpoint.  Finally,  the  air 
mass  flow  change  is  seen  to  be  small,  but  it  is  large  enough  to  produce  the 
changes  in  Tcool. 

Overall,  the  variations  in  secondary  subsystem  control  variables  are 
small  and,  except  for  one  questionable  case,  the  variables  return  to  the 
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setpoint.  The  degradation  of  response  does  not  seem  to  be  a  valid  measure  for 
the  zone  subsystem  since  its  isolated  response  was  not  designed  by  selection 
of  controller  gain  settings.  Although  the  interconnections  are  designed  to 
assist  in  maintaining  constant  zone  temperatures,  the  comparison  of  isolated 
and  primary  zone  subsystem  final  temperature  values  indicates  that  the  inter¬ 
connections  do  little  to  maintain  the  zone  temperature  when  zone  load  changes 
occur. 

In  summary,  this  chapter  has  demonstrated,  by  computer  simulation,  that 
effects  of  interconnection  of  the  isoLated  subsystem  are  of  insignificant 
magnitude.  There  are  degradations  from  the  design  responses,  but  these,  at 
worst,  increase  the  settling  time,  e.g.,  the  cooling  coil.  There  are  varia¬ 
tions  in  the  secondary  control  variables,  but  all  but  one  of  these  returned  to 
the  setpoint.  So  by  these  two  measures,  the  SISO-designed  MIMO  control  system 
is  rather  insensitive  to  the  interconnections  which  were  neglected  during  the 
design. 

It  seems  important  to  mention  here  what  was  not  seen.  The  interconnec¬ 
tion  of  the  subsystem  results  in  a  closed  loop  among  these  subsystems.  By 
examining  Figure  2,  this  loop  can  be  seen;  it  includes  the  mixing  box,  the 
cooling  coil  (or  just  the  supply  air  duct  if  the  coil  is  off),  the  zones,  and 
back  to  the  mixing  box.  Now,  when  the  return  air  fraction  is  large,  there  is 
a  continuous  loop  of  air  among  these  subsystems  further  complicated  by  the 
transport  lags  of  the  ductwork.  Problems  can  occur  if  the  zone  heat  capa¬ 
cities  are  small,  e.g.,  C  •  =  50  kJ/°C  (the  size  of  the  USA-CERL  HVAC  test 
facility  zones),  rather  than  the  size  considered  for  this  study,  Cz^  =  700 
kJ/°C  (the  average  size  of  a  zone  in  which  this  HVAC  system  might  be  used). 
Then,  this  loop  is  capable  of  sustaining  oscillations  which  do  not  decay. 

This  was  in  fact  observed  in  the  earlier  simulations  until  the  zone  size  was 
changed  from  the  smaller  to  the  larger  values  indicated.  As  long  as  larger 
zones  are  used,  this  should  not  be  a  problem. 
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VIII.  CONCLUSIONS 


The  conclusions  to  this  report  fall  into  three  categories:  modeling, 
comparison  sensitivity,  and  simulation. 

The  results  of  Chapters  III  and  IV  indicate  that  modeling  of  an  HVAC 
system  can  be  accomplished.  Prom  experimental  results,  many  of  the  responses 
can  be  described  by  single-pole  transfer  functions.  The  problem  of  hysteresis 
can  be  handled  by  the  ACSL  package.  And  finally,  the  system  pLant,  in  linear¬ 
ized  form,  can  be  described  by  a  transfer  function  matrix. 

The  comparison  sensitivity  of  the  system  could  not  be  established  by  the 
singular  value  technique  which  was  employed.  This  does  not  mean  that  the 
system  does  not  satisfy  the  comparison  sensitivity  condition,  but  that,  in 
this  instance,  the  results  of  this  test  were  inconclusive. 

The  simulations  of  the  isolated  subsystems  and  the  entire  system  were 
useful  in  establishing  two  results.  First,  the  Ziegler-Nichols  technique  for 
controller  parameter  value  selection  is  useful  for  the  PI  controllers  as  used 
on  the  three  subsystems.  But,  to  achieve  the  desired  1 /4-wave  responses,  the 
revised  Z-N  formulas  must  be  used.  Second,  by  comparing  the  isolated  subsys¬ 
tem  setpoint  step  change  responses  to  the  primary  subsystem  setpoint  step 
change  responses,  and  observing  the  associated  secondary  subsystem  control 
variable  changes,  it  was  possible  to  establish  that  the  sensitivity  of  the 
SISO  design  to  the  presence  of  the  MIMO  interconnections  was  insignificant,  as 
measured  by  the  degradation  of  the  primary  subsystem  response  and  the  varia¬ 
tion  of  the  secondary  subsystem  control  variables. 


107 


IX.  RECOMMENDATIONS 


The  recommendations  of  this  study  are  of  two  types:  I)  those  pertaining 
to  testing,  based  on  the  results  of  this  study;  and  2)  those  pertaining  to  the 
course  of  future  study. 

As  to  testing,  based  on  the  conclusions  of  this  report,  it  is  suggested 
that  if  this  design  is  also  considered  desirable  from  energy  considerations, 
it  should  be  tested  on  a  full-size  facility.  This  should  be  done  using  the 
Z-N  technique  for  controller  parameter  value  selection,  but  with  the  modified 
formulas.  Caution  should  be  taken  to  assure  that  the  zones  are  large  enough 
to  avoid  the  system  loop  oscillations  mentioned  at  the  end  of  Chapter  VII. 
Finally,  the  system  should  of  course  be  observed  during  all  seasons  before  any 
conclusions  are  made. 

For  any  future  study  utilizing  the  results  of  this  report,  two  items 
should  be  mentioned.  First,  the  model  of  the  system  might  be  improved  by 
taking  energy  considerations  into  account  when  modeling  the  fan  system, 
combining  the  E/P  and  actuator  transfer  functions,  and  finding  a  different 
description  of  the  coils.  Second,  use  the  system  plant  transfer  function 
matrix,  P' ,  to  find  feedback  and  feedforward  compensators,  G  and  H,  which 
insure  MIMO  system  criteria.  This  procedure  could  take  the  form  of  using  the 
singular  value  test  as  a  design  criterion  to  insure  comparison  sensitivity  of 
the  system. 
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